Short Range Ising Spin Glasses: a critical exponent study 
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Abstract 

The critical properties of short-range Ising spin-glass models, defined on a 
diamond hierarchical lattice of graph fractal dimension df = 2.58, 3, and 4, 
and scaling factor 2 are studied via a method based on the Migdal-Kadanoff 
renormalization-group scheme. The order parameter critical exponent (3 is 
directly estimated from the data of the local Edwards- Anderson (EA) order 
parameter, obtained through an exact recursion procedure. The scaling of 
the EA order parameter, leading to estimates of the v exponent of the cor- 
relation length is also performed. Four distinct initial distributions of the 
quenched coupling constants (Gaussian, bimodal, uniform and exponen- 
tial) are considered. Deviations from a universal behaviour are observed 
and analysed in the framework of the renormalized flow in a two dimen- 
sional appropriate parameter space. 
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In this work, we report direct estimations of the critical temperature, the order param- 
eter (3 critical exponent and the correlation length critical exponent for the short range 
Ising spin glass models on the diamond hierarchical lattices (DHL) with graph fractal 
dimension df = 2.58, 3.0 and 4.0. By direct we mean that the exponents were computed 
from the numerical values of the local Edwards- Anderson (EA) order parameter obtained, 
through an exact recursive procedure developed by the authors to study the multifractal 
structure of the local order parameter |l] . This approach is based on the Migdal-Kadanoff 
renormalization group approximation for finite hypercubic lattices @, which was proven 
to be exact, for pure systems, on the DHL. Within this method one is able to calculate 
the local magnetization of each site for each sample of the corresponding DHL as a func- 
tion of the temperature. An scaling of the order parameter at T = T c , makes it possible 
to obtain values for the exponent <fi — f3/v, and therefore the correlation-length critical 
exponent v as well. 

The spin-glass Ising model on a diamond hierarchical lattice of arbitrary graph fractal 
dimension df, is defined by the hamiltonian 

H = - JijWj (!) 

<i,j> 

where {Jij} are the random quenched nearest-neighbor coupling constants following a 
given distribution and the ex's represent the Ising variables assigned to the lattice sites. 
We consider four different distributions of coupling constants as an initial model condition, 
namely, Gaussian, bimodal, uniform and exponential, for lattices up to sixteen generations. 
The values of the Edwards- Anderson order parameter qf , defined for each lattice site, 
can be numerically evaluated for each initial distribution by means of an exact recursion 
procedure as a function of the temperature, for a general DHL with N generations fl[]. 
Adding such values over an appropriated set of lattice sites for a given sample (one 
realization of the initial set of coupling constants) and taking the average over many 
samples, we end up with the usual EA order parameter per spin, 

(2) 

j 

In the equation above, stands for the configurational average, whereas the sum 
is restricted to those lattices sites within a given shortest path joining the two end sites 
of the hierarchical lattice 0; for a DHL with N generations and scaling factor b, one 
has M = 6+1. Besides saving computer memory, this choice is expected to yield 
a stochastic representative set of local order parameters; after the average over many 
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different samples is performed, any peculiarities of a particular path become irrelevant. 
The critical exponent (3 is defined by 



q* A ~ [(T c - T)/T c f (T — > T~), (3) 

where T c is the spin-glass critical temperature, estimated by following numerically the 
probability distribution P(ty) associated with the thermal transmissivities, {ty} = 

tanh(/?{ Jy}). For a centered initial distribution = («/#) = 0), T c is defined as 

/ 9 V 1 / 2 

the temperature at which the width (tfj) remains constant under the renormalization 
process . Moreover, the correlation- length critical exponent v can be also estimated by 
scaling the order parameter at the critical point, i.e., 



q 



EA 



M /?/„ (T = T C ). (4) 



Actually, we estimate the (3 critical exponents for a dj — 2.58, 3.0 and 4.0 DHL with 
scaling factor b = 2 by considering lattices with N = 15 generations and 300 samples. 
To obtain the ratio (3/u we consider N varying from 8 to 15 and the number of samples 
varying from 1000 ( for smaller systems) to 300 (larger systems). The values of the local 
EA order parameters were calculated following the method developed by the authors fl|], 
considering four distinct initial distributions of coupling constants namely, the Gaussian, 
bimodal, exponential and uniform distributions, all presenting an unitary width, given 
respectively by 

P(J M ) = ^exp(-iJ^.) (5) 

P(Ju) = k[8{Ju-l)+8{Jij + l)) 
P(J M ) = ^exp(- v / 2|J M |) 



P(Ji 



i,3 1 



^ if ~VS < J itj < y/3 
(otherwise) 



In Table I, we display our results for (3 and v critical exponents and critical tem- 
peratures for lattices with three graph fractal dimension df. It should be mentioned that 
as far as the hierarchical lattice is concerned, the only approximations involved in the 
present approach are the finite sizes, and the number of samples investigated. 

For all lattice dimensionalities the critical temperatures were obtained by following 
the flow of the renormalized distribution width. In all cases the critical temperatures 
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are smaller as higher is the kurtosis of the initial distribution of coupling constants. For 
lattices with graph fractal dimension df = 2.58, which correspond to a DHL basic unit 
with three connections and scaling factor 2, the critical temperatures are close but different 
from zero, indicating that the SG lower critical dimension di is still lower than df = 2.58. 
Actually, within the MKRG approach de was estimated to be around 2,43 [f7ULl|. 

The critical exponents (3 and v were obtained for all considered distributions except 
of the exponential one due to an avoided numerical overflow. Nevertheless, the figures 
obtained so far, shown an universal behavior within the error bars. 

For lattices with df = 3 and 4, we have recovered recent results for the values of the 
critical temperatures obtained within the same approach. Moreover, whenever possible to 
compare, the agreement of our results with those obtained from other methods, developed 
for the cubic lattice (e.g., extensive Monte Carlo simulations high temperature series 
expansions ||) is very impressive. However, as far as the MKRG scheme is concerned 
there is a strong disagreement, in the Gaussian case, between our estimates of the expo- 
nent v (1.8), and that obtained by Southern and Young (2.78) ||, evaluated also with a 
similar unit cell. The basic differences between their approach and ours are: (i) they force 
the bond distribution into a Gaussian at each step; (ii) they compute the v exponent by 
monitoring the width of the distribution in the vicinity of the critical point, looking at 
the fluctuations around its fixed value, linearizing its recursion relation (linear approxi- 
mation). It is a well-known fact |3],|7j that the fixed-point distribution is indeed very close 
to the Gaussian one, so approximation (i) should not be too drastic. However the linear 
approximation has been criticized, and it may not be appropriate for the spin glasses, as 
it is for pure systems |§; we suspect that this is the main reason for the disagreement 
between our estimate and the one of Ref. ||. It is worth to mention that, for the case of 
a Gaussian distribution, the v exponent calculated in the present work recovers the one 
obtained by McMillan |§, applying a method, called domain- wall renormalization group, 
based on the calculation of the domain-wall free energies for very small finite-size cubic 
lattices. 

For df = 3, our results show that slight, but undeniable, different values of the ex- 
ponent P were obtained for the corresponding initial distributions of coupling constants, 



suggesting absence of universality in this model, as claimed by other authors [10|. How- 
ever, this apparent lack of universality obtained within the present MKRG scheme should 
be interpreted carefully, taking into account the nature of the flow associated with the 
evolution of the coupling constants distributions, under the renormalization process. This 
is represented schematically in Fig. 1. If the renormalized distribution is forced to be 
shape-invariant at each renormalization step (as done in Ref. H] for the Gaussian case), 
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then the renormalized flow is restricted to a one- dimensional parameter subspace. How- 
ever, if the distribution is let to evolve freely under the renormalization process, one may 
represent its flow in a higher dimensional parameter space. An illustrative two dimen- 



sion projection of this space is the plane (t^-) 1//2 versus k B T / (Jfj) 1 ^ 2 Within this 



representation, the parameter space is separated in two regions, each one governed by 
its respective attractor (spin-glass and paramagnetic stable fixed points). In the frontier 
of these regions, one finds the saddle-point spin-glass critical point, characterized by a 
fixed-point distribution, which should remain unchangeable under renormalization. Since 
the analytical form of this distribution in unknown, the genuine critical properties asso- 
ciated with the spin-glass phase transition can not be determined through the present 
approach. Therefore, the results herein presented for each distribution, should be consid- 
ered as "pseudo-critical", rather than critical properties. In Fig. 1 we also show typical 
flows starting from points slightly above and below the corresponding critical temperature 
of the bimodal distribution. Here, again, we notice that the position of the line (tjj\ 
corresponding to each considered distribution also follows the decreasing of the kurtosis 
as we go from the left to the right in Fig. 1 |f^,p|. In this work, the pseudo-critical (3 
exponents were obtained directly from the values of the local EA order parameters, calcu- 
lated in the vicinity and below the critical temperatures associated with each considered 
probability distribution. For the present DHL, where | of the total number of sites belong 
to the last generation, and whose local order parameters are calculated with the coupling 
constants introduced by the not-yet-renormalized (initial) distribution, different values of 
(3 appear, describing the behavior of the system near these critical points; this should not 
be interpreted as a breakdown of universality. To obtain the genuine critical exponents, 
one should probe the critical region around the saddle-point spin-glass critical point a task 
which deserves a careful numerical analysis. We notice that for higher space dimensions 
(i.e., large number of parallel branches in the basic unit cell), the genuine fixed-point 
distribution should approach the Gaussian one [jP3^. 

In higher dimensions (df = 4) the f3 exponent approaches to the mean field value 
for spin glasses with slight differences arising from the consideration of distinct initial 
distribution of couplings as above discussed for the three dimensional case. 

Finally, we have also estimated the v exponent by calculating the scaling exponent <p = 
P/u at the critical temperature corresponding to each distribution. We have obtained the 
same value for the correlation- length exponent, within the error bars, for the considered 
distributions. Moreover, these values are in good agreement with the corresponding values 
obtained by high-temperature series expansions and numerical simulations for the spin- 
glass Ising model on a three dimension cubic lattice . It is amazing how the 
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estimates for the spin-glass critical temperatures and exponents, on a diamond hierarchical 
lattice with scaling factor b = 2 and graph fractal dimension df — 3, can be so close to 
those of the corresponding Bravais lattice; why this happens for complex systems like spin 
glasses, and not to simple ferromagnets, is a point which is now deserving further studies 



11]. For lower and higher dimensions both exponents are slightly independent from the 



initial distribution do not suggesting breaking of universality. 
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TABLES 
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TABLE I. The estimated values for the critical temperature T c and critical exponents (3 ,/3/v 
and v . The columns are distributed accordingly with the increasing of the initial distribution 
kurtosis. Note that the critical temperatures decreases as the kurtosis increases. 
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FIGURES 




FIG. 1. Schematic flow diagram in the appropriated parameter space for the Ising SG model 
on the DHL. The full thick line represents the "fixed-point distribution", whereas the dashed 
lines, the non-renormalized (initial) distributions: exponential, Gaussian, uniform and bimodal. 
Full lines indicate the flow starting from a hypothetical point, close to critical temperature of 
the bimodal initial distribution. The full circle represents the spin-glass fixed critical point and 
open circles, the critical points corresponding to each of the initial distributions 
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